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Abstract 

CP violation in B — > J/ ipK* decays is studied using an angular analysis in a data sample of 253 
fb _1 recorded with the Belle detector at the KEKB e + e~ collider. The flavor separated measure- 
ments of the decay amplitudes indicate no evidence for direct CP violation. T-odd CP violation is 
studied using the asymmetries in triple product correlations, and the results are consistent with the 
Standard Model null predictions. The time-dependent angular analysis gives the following values 
of CP-violating parameters: sin2</>i = 0.24 ± 0.31 ± 0.05 and cos2^i = 0.56 ± 0.79 ± 0.11. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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An angular analysis of B meson decay to two vector mesons is a sensitive probe of new 
physics. There are three classes of parameters obtained through the angular analysis. The 
first is the measurement of the decay amplitudes of the three helicity states. These can be 
obtained by the time-integrated angular analysis of flavor specific decays. The comparison 
of the amplitudes between flavor separated samples probes direct CP violation. The second 
is the triple product correlation, which can be extracted from the measured decay ampli- 
tudes. This quantity is sensitive to T-odd CP violation. The third class is comprised of the 
CP parameters (sin20i and cos20i) that are measured through a time-dependent angular 
analysis. In particular, the measurement of cos20i, which appears in the time- dependent 
interference terms, is important both to solve the two-fold ambiguity in 20i and to test the 
consistency of this determination with the more precise value from other b — > ccs decays. In 
this paper, we report the measurements of these parameters for B — > J '/ ipK* decays. 

The data sample used in this analysis corresponds to an integrated luminosity of 253 
fb _1 recorded with the Belle detector jlj at the KEKB electron-positron collider 0. 
Four decay modes are reconstructed: B° — ► J/if)K*°; K*° — ► K + 7T~ and K^n , and 
B + — > J/ipK* + ; K* + -h> K + tt° and Kgir + . The charge conjugate modes are included every- 
where unless otherwise specified. The reconstruction is done using the criteria described in 
Ref. 0. A J /if) candidate is reconstructed from two oppositely charged tracks identified as 
leptons. A K* candidate is selected if the absolute difference between the invariant mass of 
an identified Kir pair and the nominal i^*(892) mass is less than 75 MeV/c 2 . Candidate B 
mesons are reconstructed by combining a J /if) candidate with a K* candidate and examin- 
ing two quantities in the center-of-mass of the T(4S): the beam-constrained mass calculated 
using the beam energy in place of the reconstructed energy (M bc ), and the energy difference 
between the B candidate and the beam energy (AE). M\> c is required to be in the range 
5.27-5.29 GeV/c 2 . For the modes with a charged (neutral) pion, the energy difference must 
satisfy \AE\ < 30MeV (-50MeV < AE < 30MeV). To eliminate slow tt° backgrounds, 
the angle 9k* of the kaon with respect to the opposite of B direction in the K* rest frame is 
required to satisfy cos9k* < 0.8. When an event contains more than one candidate passing 
the above requirements, the best combination is selected based on a \ 2 calculated using 
Mb c and AE. Figure Q shows the M\ )C distributions for the candidates of two neutral B 
decay modes within the AE signal window. After all selections are applied, the remaining 
numbers of events in the signal region of the M^ C -AE plane are 8194 for J/if>K*°(K + 7r~), 
363 for J/ipK*°(K° s ir°), 2222 for J /if)K* + (K + ^) and 2168 for J/ipK* + (K^7T+). 

The angular distribution of the products of B — > J/ipK* decays is described using three 
angles in the transversity basis |4j. The direction of motion of the J /if) in the rest frame 
of the B candidate is defined to be the x-axis. The y-axis is chosen along the direction of 
the projection of the K momentum into the plane perpendicular to the x-axis in the B rest 
frame. The z-axis is then perpendicular to the x-y plane according to the right-hand rule. 
The angle between the positive lepton (Z + ) and the z-axis in the J /if) rest frame is defined 
as 9 tr - The angle between the x-axis and the projection of the l + momentum onto the x-y 
plane is defined as 4> tr m the same frame. The angle 6k* is defined earlier. 

The distribution of the angles as a function of the decay time difference between B and 
anti-S mesons (At) is described as follows j5j: 



1 rf 4 r(fltr,&nflAr.,At) 9 e-l A *l/ r * 

T d cos 6 tr d(f>trd cos 9k* dAt 32n 2tb 
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FIG. 1: M bc distributions for the candidates of (a) B° -> J/if>K*° (K + ir~), and (b) B° -> 
J/tpK*°(KgTr°). The solid line shows the projection of the two-dimensional fit to M]-, c and Ai£ 
while the dashed line is the estimated background contamination. 



where the angular terms gi are denned as 

gi = 2 cos 2 9k* (1 — sin 2 9 tr cos 2 tr ) 

g% = sin 2 9k* (1 — sin 2 tr sin 2 tr ) 

gs = sin 2 sin 2 9t r 

g± = — (l/y/2) sin 29k* sin 2 ^ r sin 20f r 

(?5 = sin 2 sin 2# tr sin <fi tr 

ge = (1/V2) sin 29 K * sin 29 tr cos <p tr , 

and the amplitude terms di as 

di = |t4 | 2 (1 + rj sin 20i sin Am At), 
a.2 = \A\\ 1 2 (1 + 77 sin 20i sin Am At), 
a 3 — I 1 2 ( 1 — sin 20i sin Am At), 
a 4 = Re(A|v4 )(l + r/sin20isin AmAt), 
05 = r]lm(A*^A^) cos AmAt 

-r]Re(A*A ± ) cos 20 a sin AmAt, 
a 6 = ^Imf^A^) cos AmAt 

— r]Re(y4Q74^) cos 20i sin AmAt. 

Here, Aq, An and A± are the complex decay amplitudes of the three helicity states in the 

transversity basis, and 77 = +1 (—1) for B° or £> + (B° or £>~). Am is the B° — 5° mixing 
parameter, which is zero for charged B meson decays, and tb is the lifetime of a B meson. 
T is the decay rate to each final state. Two CP violation parameters appear in the formula, 
viz, sin 20! and cos20i. They can take non-zero values for the decay into CP eigenstate 
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B° — > J / ip K*° (K s tt°) that occurs through the same quark level diagram as that of the 
golden mode for sin20i measurement, B° — > J/ipK s . 

The determination of the decay amplitudes and CP parameters are performed using an 
unbinned maximum likelihood method, taking into account the detection efficiency and 
backgrounds. The probability density function (PDF) for an event is defined as 

V = fsig{M hc ,AE)e(6 tr ,(f) tr ,8 K *)- 



r d cos 9 tr d(f)trd cos 9x*dAt 

e -\M/r B f . 
+ ^ < J2 /c/( M bc, AE)A cf (6 trj fa, 9 K *) 

+ fnr(M\y C , AE)A nr (0tr, fa, ®K* 

+6(At)f cb (M hc , AE)A cb (9 tr , fa, 6 K .), (2) 

where f s i g , f c f, f nr and f cb are the respective fractions of signal, cross- feeds, non-resonant 
production, and combinatorial background components as functions of AE and M bc , while 
A c f, A nr and Ado are the corresponding angular shape functions, and e is the detection 
efficiency function. 

Three separate background sources are considered in the fit: cross- feeds that are contam- 
inations from other K* subdecays, non-resonant production of Kir including contaminations 
from higher resonance tails, and combinatorial background. The fraction and angular shape 
of cross-feeds are estimated from Monte Carlo. The fraction of non-resonant production 
is estimated to be 6.8% from a fit to the invariant mass distribution of Kir pairs in an 
alternate data sample that is assembled without the K* selection criteria. The fit is per- 
formed with two Breit-Wigner functions for the X*(892) and if|(1430) resonances and a 
threshold function describing the non-resonant production, taking into account the con- 
taminations of cross-feeds and combinatorial background. The phase space factor includes 
a J/ip recoil correction. The angular shape is obtained from samples in the mass region 
1.0GeV/c 2 < M(Ku) < 1.3GeV/c 2 . 

The fractions of signal and combinatorial background are estimated from a fit to the 
samples in the region of 5.2 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.1 GeV < AE < 0.1 GeV. 
The angular shape of combinatorial background is obtained from the sub-sample with 
5.2GeV/c 2 < Mbc < 5.26GeV/c 2 used in the fit above. The fractions of both signal 
and backgrounds are parameterized as two-dimensional functions of M^ c and AE. The 
background angular shapes are parameterized for each of three angles separately. For each 
signal mode, the detection efficiency is parameterized as a three-dimensional function whose 
parameters are obtained by a fit to a high-statistics Monte Carlo sample. The function is 
almost flat except in the region cos 9k* ~ 1, where the pion is slow so that the efficiency is 
reduced. 

The decay amplitudes are determined by fitting the time-integrated angular distribution 
to three measured angles. Eqs.^andElare integrated over At, where terms with 6 and 
cos Am At become unity while terms with sin Am At become 0. The value of r\ is determined 
from the charge of the kaon for K* decays with a K + or of the pion with a Kg. The decay 
mode B° — ► J/tljK*°(K s 7t°) is not used, since 7] cannot be determined by this prescription. 
In the fit, the imaginary part of Aq is defined to be zero since the overall phase of the 
decay amplitudes is arbitrary. The values of the other five parameters, \A \ 2 , \A\\\ 2 , \A^\ 2 , 
arg(v4||) and arg(Aj_), are determined in the fit. There is a two fold ambiguity in arg(Aii) and 
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TABLE I: Measured decay amplitudes for B° and B 
the second is systematic. 



decays. The first error is statistical while 
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B u 



B 



B" 



B 



B+ 



B~ 



B~* 



B- 



0.571 ±0.015 
0.216 ±0.017 
0.213 ±0.017 
-2.934 ±0.134 
2.878 ± 0.088 



0.578 ±0.016 
0.244 ±0.018 
0.178 ±0.017 
-2.851 ± 0.114 
2.993 ± 0.089 



0.574 ±0.012 ±0.009 
0.231 ±0.012 ±0.008 
0.195 ±0.012 ±0.008 
-2.887 ±0.090 ±0.008 
2.938 ±0.064 ±0.010 



0.600 ±0.020 
0.194 ±0.019 
0.206 ±0.019 
-3.070 ±0.142 
2.964 ± 0.099 



0.608 ±0.021 
0.243 ±0.021 
0.149 ±0.019 
-3.129 ± 0.172 
2.988 ±0.121 



0.604 ±0.015 ±0.018 
0.216 ±0.014 ±0.013 
0.180 ±0.014 ±0.010 
-3.090 ±0.108 ±0.006 
2.983 ± 0.076 ± 0.004 
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FIG. 2: Distributions of projected angles for B° — > J/ipK* (K + ir~). Solid lines show results of the 
fit. The data points are corrected for the detector efficiency and the backgrounds are subtracted. 

arg(A^) We take the choice consistent with the s-quark helicity conservation hypothesis, 
which is shown by BaBar to be the physical choice |7[. The normalization condition of the 
amplitudes, \A \ 2 + \A\\\ 2 + \A±\ 2 = 1, is taken into account by adopting the extended 

likelihood defined as — InL = — X^i s 111 Qi + -^exp — hi(iV ex p), where Qi is the value of 
the PDF for each event, and iV b s is the number of events used for the fit. iV ex p is defined to 
be A' r b s -(|^4o| 2 +|^4||| 2 + |^4±| 2 ) to incorporate the normalization condition. The normalization 
of the PDF is recalculated whenever the fit parameters change. The two charged B decay 
modes are combined by defining a single likelihood. 

The decay amplitudes determined from the fit are summarized in Table HI The obtained 
values are consistent between the two flavors both in neutral and charged B decays, indicat- 
ing no evidence for direct CP violation. The flavor averaged values are consistent with our 
previous measurement and that by BaBar using 83 fb _1 0]. Small discrepancies from tt 
are observed in arg(Ay) and arg(y4jj for both B° and B + decays. The difference of these 
two phases is 0.458 ± 0.110 rad in B° decays, which is shifted from by more than 4a. 
This is interpreted as evidence for the existence of final state interactions. Fig. El shows the 
projected angular distributions for B° — > J/i/jK* (K + it~) decays. 

Systematic uncertainties in the fit are determined for: 1) detection efficiency (MC statis- 
tics and effect of polarization), 2) background angular distribution functions, 3) background 
fractions, 4) slow pion efficiency, and 5) non-resonant decay polarization effect. The effect 
of the uncertainty in the fraction of the non-resonant production is estimated by varying 
the value by ±5cr to take into account the possible contamination of other resonance tails. 
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The uncertainty in its angular shape is determined by comparing with results that assume 
phase space decay. 

The triple product correlations in B — > J/ipK* decays take the form of p ■ (v\ x v 2 ) in 
the B rest frame, where p is the momentum of J /if} or K*, and Vi(v 2 ) is the polarization 
vector of J/if>(K*) [8]. They are odd under time reversal (T-odd) and their asymmetries 
are sensitive to direct CP violation even when the strong phase difference is small. The 
asymmetries are defined using the decay amplitudes as: 

(1) _ lm(A ± A*) (2) _ MA ± Al) 

a t Al + Aj + AY Al + A\ + A\- [> 

The corresponding asymmetries for anti-i? decays are defined as A^ and . The Standard 
Model predicts tiny values for these asymmetries and no difference between B and anti- 
B mesons. Substituting the measured amplitudes in Eq. El the obtained triple product 
asymmetries are listed in Table |H] As seen, all the obtained asymmetries are small. Tiny 
discrepancies from zero are considered to be due to final state interactions. No difference 
between A^ and Ay nor Ay and Ay is observed, which is consistent with the absence of 
T-odd CP violation. 

TABLE II: Measured asymmetries in the triple product correlations. The first error is statistical 
while the second is systematic. 





B° 


B+ 


AM 


0.091 ± 0.034 ± 0.007 


0.062 ±0.038 ±0.005 




-0.098 ±0.032 ±0.003 


-0.049 ±0.034 ±0.002 


J(l) 


0.047 ± 0.031 ± 0.007 


0.046 ± 0.039 ± 0.005 


1(2) 


-0.089 ±0.029 ±0.003 


-0.031 ± 0.039 ± 0.002 


w^m 


0.044 ± 0.046 


0.016 ±0.054 


\ A ® 1(2), 


0.009 ± 0.043 


0.018 ±0.052 



The CP violation parameters sin 20! and cos 2<pi are determined by fitting the time- 
dependent angular distribution in Eq. ^ to the measured angles and At simultaneously 
in the CP decay mode B° — > J/if>K*°(Kg7i°) using the PDF in Eq. |21 The procedures to 
measure At and to determine the flavor of the decaying B° meson are described elsewhere (§[ . 
The flavor tagging procedure gives the flavor q of the tag-side B meson, where q = ±1(— 1) 

for B°(B ), and the probability w that this flavor determination is incorrect. The value of 
r\ is given by —q(l — 2w) within the signal angular distribution in Eq. ^ The difference in 

w between B° and B mesons is also considered. 

Each term in Eq. |21 is convolved with the appropriate resolution functions separately 
for the signal, backgrounds having the B° lifetime (namely, cross-feeds and non-resonant 
production), and the combinatorial background with a zero-lifetime 5-function shape. The 
resolution functions are obtained from fits to the At distributions measured for various 
data samples. Null CP asymmetry is assumed for the backgrounds. In the fit, the decay 
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FIG. 3: Raw asymmetry in measured At between samples tagged as q = +1 and q = — 1. The 
solid line shows the projection of the fit. 



amplitudes are fixed at the values obtained for B° decays. The lifetime and mixing parameter 
are set to PDG values jioj . From the fit to the data, we obtain 

sin 20i = 0.24 ±0.31 ±0.05, 
cos20! = 0.56 ±0.79 ±0.11. 

When we fix the value of sin20i to the world average value (0.726) [IH, the value of cos20i 
becomes 0.87 ± 0.74 ± 0.12. The positive sign of cos20i is consistent with the measurement 
by BaBar 0; however, we cannot exclude negative values with our current statistical errors. 
The raw asymmetry in the measured At distribution between samples with q = ±1 and —1 
is shown in Fig. El with the projected result of the fit. 

Systematic uncertainties in the fit are determined in the same manner as those in the 
b — > ccs sin 20! measurement 0. In addition, the uncertainties that come from the angular 
analysis are estimated similarly as that in the decay amplitude measurement. The possible 
bias in the fit is checked by applying the same fitting procedure to the sample of B — > 
J/ipK*°(K + 7c-) decays. We obtain "sin20i" = -0.047 ± 0.067 and "cos20i" = -0.111 ± 
0.161, which are consistent with zero as expected. 

In summary, a full angular analysis is performed for B — > J/ipK* decays. The complex 
decay amplitudes are measured by a simultaneous fit to three transversity angles. The 
measured values are consistent between the two B flavors both in neutral and charged B 
meson decays, and no direct CP- violating effect is observed. The difference of arg(A||) and 
arg(AjJ of B° decays is shifted from by more than 4a, which is interpreted as evidence 
for the existence of final state interactions. The differences between the asymmetries of 
triple product correlations for B and anti-5 mesons are consistent with zero in both neutral 
and charged B mesons, and no T-odd CP-violating new physics effect is observed. The 
time-dependent angular analysis performed for B° — > J/tpK*° (K*° —>■ K^tt ) decays gives 
the CP violation parameters sin20i = 0.24 ± 0.31 ± 0.05 and cos20i = 0.56 ± 0.79 ± 0.11. 
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Fixing sin 20i at the world average value (0.726) gives cos20i = 0.87 ±0.74 ±0.12. The sign 
of cos 20i is positive, although we cannot exclude negative values with current statistical 
errors. 
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